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Abstract 

We investigate to which extend a THz laser pulse can be used to pro- 
duce field-free molecular orientation at high temperature. We consider 
laser pulses that can be implemented with the state of the art technology 
and we show that the efficiency of the control scheme crucially depends 
on the parameters of the molecule. We analyze the temperature effects 
on molecular dynamics and we demonstrate that, for some molecules, a 
noticeable orientation can be achieved at high temperature. 

1 Introduction 

Manipulating the molecular rotational degree of freedom remains a goal of pri- 
mary interest in photochemistry with applications extending from chemical re- 
activity to nanoscale design yUi2ii3j. In this framework, molecular ahgnment and 
orientation constitute a well-established topic both from the experimental and 
theoretical points of view [H [5] . While the alignment process is now well under- 
stood in the adiabatic [B] or sudden regime |7J with recent extensions such as the 
molecular classical rotation [S], the deflection of aligned molecules [5], the pla- 
nar alignment |10l [TTl \T2\ or the analysis of the dissipation effects [13' 13] , work 
remains to be done in order to control and produce molecular orientation with a 
sufhcient high efficiency. On the theoretical side, several basic mechanisms have 
been proposed, built on intuitive or optimal control strategies [TB] . Among 
others, we can cite the kick mechanism which consists in a sudden impact to 
the molecule by a half-cycle pulse (HCP) [13 [H [H [501 [HI 122 , its combination 
with a laser field [Ml HI PS] or the {uj~2lj) scheme [?7l 1^ 1^ [5^1 [gT| [521 . 
Due to the efficiency of the first process based on its asymmetric temporal shape, 
most of the theoretical works have focused on its application [341 [551 [551 [571 [55] . 
However, recent studies [391140] have pointed out the inherent experimental diffi- 
culties associated to the use of such pulses, which are distorted when they prop- 
agate in free space as well as through focusing optics [41J. This phenomenon, 
due to the DC component of the field, makes thus problematic the experimen- 
tal implementation of these techniques in the control of molecular rotation. In 
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this framework, a fundamental question is whether it is possible to orient linear 
molecules in the THz regime by using only zero area laser pulses. These latter 
do not contain DC field and are therefore free of these propagation distortions. 
This problem has been recently addressed theoretically [IDl H5] and 

experimentally [15], but no systematic study of the efficiency of this process 
has been done. Note that this question is not trivial since the sudden impact 
approximation predicts no post- pulse orientation in this regime [T^f^UlHO] . We 
present in this paper a complete analysis of this control strategy by consider- 
ing different linear molecules. We establish under which conditions this process 
is efficient and we analyze its robustness with respect to temperature effects. 
Two mechanisms leading to molecular orientation are identified. The first one 
is valid at low temperature, while the second process is only efficient for higher 
temperatures. In this second non-intuitive control scheme, we show the positive 
role of temperature effects in the orientation mechanism. 

The paper is organized as follows. The model system is presented in Sec. 
[2l The different numerical results are discussed in Sec. [3l Conclusions and 
prospective views are given in Sec. H) 



2 The model system 

We consider the control of a linear polar molecule by a linearly polarized THz 
laser field E{t) of zero area. The molecule is assumed to be in its ground vibronic 
state. Within the rigid rotor approximation, the Hamiltonian of the system can 
be written as 

H{t) = BJ^ - ^ioE{t) cos 0, (1) 

where B is the rotational constant, the angular momentum operator and /io 
the permanent dipolar moment. We neglect in this paper the effect of the polar- 
izability components since the maximum intensity of the electric field remains 
moderate. The units used throughout the paper are atomic units unless other- 
wise specified. At non zero temperature, the system can be either described by 
a density matrix p(t) with a dynamics governed by the von Neumann equation 

^^ = lH{t),pit)], (2) 

where p{0) is the canonical density operator at thermal equilibrium, or by a set 
of wave functions lipJo.Moit)) satisfying each the Schrodinger equation: 

« gl = H(t)\lPj„^Mo), (3) 

with as initial state \ipjg,Mo{t = 0)) = | Jq, Mq). The second representation given 
by Eq. ([3]), which is more suited to the control mechanisms interpretation, will 
be used in this work. The expectation value (cos 6) defined by 

^ +OC A/o = Jo 

{cos d) (t) = - Cj„ (V'Jo,Mo|cOs6'|V'Jo,A./o), (4) 

Jo=0 Mo = ~Jo 



is taken as a quantitative measure of orientation, with the weights cj = e Bj{j+i)/{kBT) 

•\oo \ J 



and the partition function Z = X^jLo X]m=-j '^•/j where T is the temperature 
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Table 1: Molecular parameters of different molecules used in the numerical 
computations. Numerical values are taken to be Eq = 2.19 MV/cm, (5 = 5 ps 
and / = 0.5 THz. 



Molecule 


OCS 


HF 


LiH 


CO 


LiCl 


B (cm-i) 


0.2029 


20.956 


7.513 


1.931 


1.345 


Ho (debye) 


0.712 


1.820 


5.88 


0.112 


6.33 


A 


117.8497 


2.9167 


26.2842 


1.9479 


158.0563 


F 


13.0823 


0.1267 


0.3533 


1.3746 


1.9735 


D 


0.1911 


19.7371 


7.0760 


1.8187 


1.2668 



and kB the Boltzman constant. In the rest of the paper, we will need to separate 
the zero temperature contribution to the thermal one, which are respectively 
denoted by (cos0)o and {cos 9)t- In the zero temperature response, the sum of 
Eq. is carried out only for Jq — 0, while for the thermal contribution, the 
expectation value is computed over the other values of Jq. 
The electric field is assumed to be of the form 

Bit) = Bo fit) - Bo cos2(7r^) sin(27r/t), t e [-6/2, 6/2], 

B{t) — otherwise, 

where Bo is the amplitude of the pulse, S its duration and / its central frequency. 
By symmetry, this field has a zero area for any values of d and /. Note that the 
sudden impact approximation can be applied if the pulse duration 6 is small with 
respect to the rotational period Tper — n/B. This also means that noticeable 
orientation can be produced only for sufficient large values of 5 [T9l |40] . 
The Schrodinger equation ([3]) can be written as 

= [J' -Acosef{r)]\ijj,^Mo), (5) 
where the new dimensionless parameters are defined by: 

r = Bt, A = i^, F = l, D = BS, f 

B ' B' ' B ' 

with /(r) = cos^(7r-^) sin(27rFr). In these coordinates, the rotational period 
becomes Tper — tt and the pulse duration is D. These effective parameters 
completely describe the dynamical evolution of the system and the degree of 
molecular orientation produced. They will give a general understanding of the 
orientation mechanism free of any particular molecule. Note that the rotational 
constant _B is a crucial parameter in these new coordinates since the effective 
field, frequency and duration (A, B, D) depend on B. For numerical applica- 
tions, we will consider the parameters listed in Table [TJ 

3 Numerical Results 

We begin our study by a general analysis of the maximum degree of orientation 
that can be reached as a function of the rotational constant B and the tempera- 
ture T . A fictive molecule with a permanent dipolar moment of 1 debye has been 
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considered in the computation. The field amplitude is assumed to be Eq ~ 2.19 
MV/cm (i.e. with a peak amplitude of 2 MV/cm), which corresponds roughly 
to the maximum amplitude of THz pulses actually available experimentally [47J . 
This dependance is shown in Fig. [T] where two zones of high orientation can be 
clearly distinguished: an upper zone, denoted (I), associated to high values of 
B and low temperature and a flat zone, denoted (II) requiring small rotational 
constants and larger temperatures up to 250 K. If the rotational constant is 
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Figure 1: (Color online) Maximum orientation as a function of the rotational 
constant B and the temperature T. The field parameters are S = 5 ps, f — 0.5 
THz and Eq = 2.19 MV/cm. Two regions (I) and (II) of the diagram are 
delimited by arcs of ellipses. 



lower than 2 cm~^, we first notice that a very weak orientation is produced at 
T = K. The rotational period varying as the inverse of B, this result can be 
explained through the sudden impact approximation [121 UHl EH] ■ Another stan- 
dard feature shown in previous studies is the detrimental effect of temperature 
on molecular orientation [SJ UHl 120! ■ This behavior can be recovered in the case 
of region (I) where almost no orientation is obtained for a temperature larger 
than 100 K. By comparison, the orientation observed in zone (II) is rather un- 
expected, since the temperature effect becomes positive, no orientation being 
produced at very low temperature. One of the goal of this work will be to ex- 
plore the basic mechanism at the origin of this non trivial phenomenon. 
The orientation mechanism. 

We first present a general spectral analysis of the control problem. We use in 
the following the dimensionless coordinates introduced in Sec. [21 for which 
the molecular spectrum is a discrete spectrum composed of the frequencies 
ujj — 2(J + 1). To simplify the discussion, we associate to each frequency a 
weight P which is defined as the average of the initial thermal population of the 
levels J and J + 1, i.e. P = {cj + cj+i)/2. The spectrum of the control field is 
proportional to the one of the function /(r) and can roughly be approximated 
by a gaussian spectrum centered aX v — F with a bandwidth proportional to 
1/D. Different spectra are represented in Fig. [21 and [31 for different values of D 
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and F, respectively. In the case of Fig. [51 only the width of the field spectrum 
is changed, while only the central frequency is modified in Fig. [31 From a qual- 
itative point of view, the overlap between the molecular and the field spectra is 
a necessary condition to produce molecular orientation. Figure [H displays two 
cases with different values of D, i.e. a pulse with a large spectrum overlapping 
many transitions including the first one, and a pulse with a larger width D 
corresponding to a narrower spectrum. The first pulse can produce orientation 
at zero and non-zero temperatures, while the second one is only efficient at 
non zero temperature. This behavior is confirmed in the top panel of Fig. [31 
Figure [31 depicts two different situations, a first one with a low main frequency 
which overlaps with the first frequency spectrum and a second case where the 
main frequency is far from the first molecular transition. The first pulse should 
work at zero and nonzero temperatures due to the overlap with other transitions 
while in the second example, we cannot expect orientation at zero temperature 
but only for T ^ 0. The corresponding orientation responses are given in the 
bottom panel of Fig. [4l Comparing in Fig. [H the efficiency of the pulses of 
Fig. [21 and [31 one observes that at T 7^ 0, a stronger orientation is obtained in 
the first case of Fig. ^ due to a larger overlap of the field spectrum with the 
population distribution. In other words, a noticeable orientation is produced 
only for a large bandwidth of the laser field. 
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Figure 2: (Color online) Plot of the spectral distribution of the field at T = 
and T = 50 (T = f43.9 K for _B = 2 cm~^) on the top and bottom panels, 
respectively. The amplitude of the field spectrum is arbitrary. Control fields 
parameters are taken to he D = 1 {S = 2.654 ps for i? = 2 cm^^) for the blue 
(black) line and D = 3 {S = 7.963 ps for B — 2 cm^^) for the green (gray) line. 
Other parameters are fixed to F = 2 (/ = 0.753 THz for S = 2 cm~^) and 
A = A. The columns represent the population distribution P = (cj + cj+i)/2 at 
the corresponding frequency oj — 2(J+ 1). The quantities P and uj are unitless. 

The interpretation of the previous results leads to two basic mechanisms 
governing the orientation response. The orientation created at T ~ K can be 
associated to a rotational ladder climbing mechanism from the ground energy 
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Figure 3: (Color online) Same as Fig. [Hbut for F = 0.5 (/i = 0.188 THz for 
B = 2 cm-i) in blue or dark and F = 2.5 (/ = 0.942 THz for B = 2 cmT^) in 
green or gray. The other parameters are given hy D — i {5 — 7.96 ps for B ^ 2 
cm~^) and A = 4. 
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Figure 4: (Color online) Time evolution of (cos 9) for the cases of Fig. [5] (top 
panel) and [3] (bottom panel). The same color code as in Fig. [5] and [3] has been 
used. Solid and dashed lines depict respectively the orientation at T = and 
T — 50. The field is on for negative times and switches off at r = 0. The 
variable t is unitless. 
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level, which consists here in the successive excitation of neighboring rotational 
levels. This process is efficient only if the field spectrum allows to excite the 
first rotational frequencies with a sufficient high intensity. Note that this con- 
trol scheme has already been identified in the literature to produce molecular 
orientation |f 5| . The rotational population distribution being displaced to high 
J levels with increasing temperatures, this mechanism then loses its efficiency 
very fast. This orientation will be called a zero-temperature orientation and can 
be quantitatively measured by the partial expectation value (cos 9)o. In the case 
of an initial thermal distribution of rotational states, a new mechanism occurs. 
This situation corresponds to molecules with small values of -B, i.e. with a quite 
large frequency F. As can be seen in Fig. [3] for F = 2.5, there is no overlap with 
the population distribution at low temperature and thus no orientation. In this 
second mechanism at high temperature, the control field excites simultaneously 
many rotational frequencies (i.e. with a large laser bandwidth), creating several 
rotational waves packets IV'Jc-A/o (0) which interfere constructively to produce 
a noticeable orientation, denoted thermal orientation. This latter can be com- 
puted through the expectation value {cos9)t- In this scenario, the orientation 
reaches a maximum at a temperature different from zero and presents a slow 
decrease with increasing temperatures. 

These different conclusions can be checked in Fig. [S]where the zero-temperature 
and the thermal orientation responses have been plotted. One clearly sees in 
this figure that the zone (I) is associated to (cos0)o, where only the initial state 
0, 0) is considered. The efficiency of THz pulses in the region (II) can be inter- 
preted as a thermal orientation such that {cos 6) ~ {cos9)t- Note that the sum 
of the two figures [5] does not give exactly the result of Fig. [T]due to destructive 
interferences between the two orientation responses. 

The transition from a regime with thermal orientation to a regime with a 
zero-temperature orientation can be understood from the definition of the pa- 
rameters {F, D). Decreasing the value of B is equivalent to increase the effective 
frequency F and to decrease the temporal width D of the pulse. Starting from 
the zone (I) where the spectrum of the field overlaps the first molecular tran- 
sition, a decrease of the rotational constant B will shift the field spectrum to 
higher frequencies F. The overlap with the first frequency is removed and no 
zero temperature orientation is possible. At the same time, as a consequence 
of the increase of the width of the field spectrum, more and more molecular 
transitions can be excited by the pulse, leading thus to a thermal orientation. 

The evolution of the orientation for the molecules of Table [T] is displayed in 
Fig. \6\ All the molecules considered belong to the zone (I), except for the CO 
and OCS molecules which present an orientation close to at T ~ OK and an 
orientation of the order of 0.15 and 0.1 at T 2± lOK and T ~ 150K, respectively. 
The efficiency of the control scheme at high temperature becomes even better 
than the one obtained for other molecules such as HF or LiH, which is very 
good only at low temperature. Note that the CO and LiCl molecules have quite 
the same parameters F and D. The difference of orientation response comes 
from the difference of two order of magnitude in their effective amplitude A. 
A ladder climbing process at T = K is possible for LiCl since the tails of 
the spectrum are sufficient, while it is not the case for CO. The two molecules 
present however a thermal orientation. We also point out that the orientation 
of OCS at high temperature is due to its low value of D (very broad spectrum) 
combined with a high value of F and a large effective amplitude of the field. 
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Figure 6: (Color online) The top panel represents the evolution of the maximum 
of orientation as a function of the temperature for the molecules of Table [T] 
The field parameters of Fig. [1] have been used. The contribution of the zero 
temperature (red or dark gray) and thermal (green of light gray) wave packets 
have been plotted for the LiCl (solid lines) and HF (dashed lines) molecules in 
the bottom panel. The blue (black) lines corresponds to the total orientation 
response. The temperature T is expressed in Kelvin. 
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In Fig. [SI note the smooth evolution of the orientation, except in the cases of 
the LiH and LiCl molecules where a slope change occurs. This feature can be 
explained by a transition from a zero-temperature orientation to a thermal one 
as shown in Fig. ^ (bottom). It can be seen that the non smooth point of the 
curve for the LiCl molecule can be viewed as a limit point for which the thermal 
orientation becomes predominant with respect to the zero-temperature one. 




0.5 1 

Logio(T) 



Figure 7: (Color online) Maximum orientation as a function of the field am- 
plitude Eq and the temperature T. The other field parameters are (5 = 5 ps 
and f — 0.5 THz. The LiCl molecule is considered for the computations. The 
temperature T is expressed in Kelvin. 

We finally explore the dependance of the orientation with respect to the 
field strength. A global description of this sensitivity as a function of the tem- 
perature is displayed in Fig. [71 At zero temperature, we observe a quadratic 
increase of the orientation up to Eq ~ 0.6 MV/cm. Higher amplitudes lead 
to a more chaotic evolution of the orientation characterized by the occurrence 
of maxima and minima. At high temperature, the same erratic distribution of 
the orientation response can be seen with, however, a larger periodicity. We 
thus conclude that the molecular orientation is less sensitive to amplitude field 
changes at non-zero than at zero temperature. We also observe that the thermal 
orientation can be produced at a lower intensity than the zero temperature one. 



4 Conclusion 

This paper has focused on the use of THz laser pulses for controlling the orien- 
tation dynamics of linear molecules. Numerical tests have shown the efficiency 
of the proposed control scheme, even at high temperature for some molecules. 
Until now, most of the works have envisaged the production of molecular ori- 
entation with short laser pulses characterized by a non zero time average of 
the electric field as, e.g., with the use of HCPs. Indeed, the sudden impact 
approximation predicts no post-pulse orientation when this average vanishes. 
Opposite to this accepted fact, we have shown that a significant orientation can 
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be obtained in two different regimes. Tfie first one corresponds to a standard 
situation with low temperature and liigh rotational constants. In the second 
case, for an adequate choice of the pulse parameters, we have established that 
the temperature plays an active role in the production of molecular orienta- 
tion. This study calls for further experimental investigation of the use of such 
laser pulses in order to complete the initial work of Ref. In particular, 

one objective could be to demonstrate experimentally the existence of thermal 
orientation in molecules such that OCS or CO. 
Acknowledgment 

We are grateful to E. Hertz for discussions. 

References 

R. C. C. Brif and H. Rabitz, New J. Phys. 12, 075008 (2010). 

M. Shapiro and P. Brumer, Principals of quantum control of molecular 
processes (Wiley, New York, 2003). 

S. Rice and M. Zhao, Optimal control of molecular dynamics (Wiley, New 
York, 2000). 

H. Stapelfeldt and T. Seideman, Rev. Mod. Phys. 75, 543 (2003). 

T. Seideman and E. Hamihon, Adv. At. Mol. Opt. Phys. 52, 289 (2006). 

B. Friedrich and D. Herschbach, Phys. Rev. Lett. 74, 4623 (1995). 

M. Leibscher, I. S. Averbukh, and H. Rabitz, Phys. Rev. Lett. 90, 213001 
(2003). 

S. Zhdanovich, A. A. Milner, C. Bloomquist, J. Flob, I. S. Averbukh, J. W. 
Hepburn, and V. Milner, Phys. Rev. Lett. 107, 243004 (2011). 

E. Gershnabel and I. S. Averbukh, Phys. Rev. Lett. 104, 153001 (2010). 

M. Z. Hoque, M. Lapert, E. Hertz, F. Billard, D. Sugny, B. Lavorel, and 
O. Faucher, Phys. Rev. A 84, 013409 (2011). 

M. Lapert, E. Hertz, S. Guerin, and D. Sugny, Phys. Rev. A 80, 051403 
(2009). 

D. Sugny and M. Joyeux, J. Chem. Phys. 112, 31 (2000). 

S. Ramakrishna and T. Seideman, Phys. Rev. Lett. 95, 113001 (2005). 

D. Sugny, C. Kontz, and H. Jauslin, Phys. Rev. A 74, 053411 (2006). 

J. Salomon, C. M. Dion, and G. Turinici, J. Chem. Phys. 123, 144310 
(2005). 

M. Lapert, R. Tehini, G. Turinici, and D. Sugny, Phys. Rev. A 78, 023408 
(2008). 

R. R. Jones, D. You, and P. H. Bucksbaum, Phys. Rev. Lett. 70, 1236 
(1993). 



11 



[18] T. Seideman, J. Chem. Phys. 115, 5965 (2001). 

[19] C. M. Dion, A. Keller, and O. Atabek, Eur. Phys. J. D 14, 249 (2001). 
[20] N. E. Henriksen, Chem. Phys. Lett. 312, 196 (1999). 
[21] M. Machholm, J. Chem. Phys. 115, 10724 (2001). 

[22] C. M. Dion, A. B. Haj-Yedder, E. Cances, A. Keller, C. L. Bris, and 
O. Atabek, Phys. Rev. A 65, 063408 (2002). 

[23] E. Gershnabel, I. S. Averbukh, and R. J. Gordon, Phys. Rev. A 74, 053414 
(2006a). 

[24] E. Gershnabel, I. S. Averbukh, and R. J. Gordon, Phys. Rev. A 73, 063402 

(2006b). 

[25] D. Daems, S. Guerin, D. Sugny, and H. R. Jaushn, Phys. Rev. Lett. 94, 
153003 (2005). 

[26] S. D. L Znakovskaya, F. A. D. Ray, N. G. Johnson, I. A. Bocharova, M. Ma- 
grakvehdze, B. D. Esry, C. L. Cocke, I. V. Litvinyuk, and M. F. Khng, Phys. 
Rev. Lett. 103, 153002 (2009). 

[27] R. Tehini and D. Sugny, Phys. Rev. A 77, 023407 (2008). 

[28] S. Zhang, C. Lu, T. Jia, Z. Wang, and Z. Sun, Phys. Rev. A 83, 043410 

(2011a). 

[29] R. Tehini, M. Z. Hoque, O. Faucher, and D. Sugny, Phys. Rev. A 85, 043423 

(2012). 

[30] S. Zhang, C. Lu, T. Jia, Z. Wang, and Z. Sun, J. Chem. Phys. 135, 034301 
(2011b). 

[31] S. Zhang, J. Shi, H. Zhang, T. Jia, Z. Wang, and Z. Sun, Phys. Rev. A 83, 
023416 (2011c). 

[32] J. Wu and H. Zeng, Phys. Rev. A 81, 053401 (2010). 

[33] Y. Hyeok, K. H. Taek, K. C. Min, N. C. Hee, and L. Jongmin, Phys. Rev. 
A 84, 065401 (2011). 

[34] I. S. Averbukh and R. Arvieu, Phys. Rev. Lett. 87, 163601 (2001). 

[35] D. Sugny, A. Keller, O. Atabek, D. Daems, C. M. Dion, S. Guerin, and 
H. R. Jauslin, Phys. Rev. A 69, 033402 (2004a). 

[36] D. Sugny, A. Keller, O. Atabek, D. Daems, C. M. Dion, S. Guerin, and 
H. R. Jauslin, Phys. Rev. A 72, 032704 (2005a). 

[37] D. Sugny, A. Keller, O. Atabek, D. Daems, C. M. Dion, S. Guerin, and 
H. R. Jauslin, Phys. Rev. A 71, 063402 (2005b). 

[38] C.-C. Shu, K.-J. Yuan, W.-H. Hu, J. Yang, and S.-L. Cong, Phys. Rev. A 
78, 055401 (2008). 



12 



[39] J. Ortigoso, arXiv p. 1108.3991 (2011). 

[40] D. Sugny, A. Keller, O. Atabek, D. Daems, S. Guerin, and H. R. Jauslin, 
Phys. Rev. A 69, 043407 (2004b). 

[41] D. You and R H. Bucksbaum, J. Opt. Soc. Am. B 14, 1651 (1997). 

[42] C.-C. Shu, K.-J. Yuan, W.-H. Hu, and S.-L. Cong, J. Chem. Phys. 132, 

244311 (2010). 

[43] C.-C. Shu, K.-J. Yuan, W.-H. Hu, and S.-L. Cong, Phys. Rev. A 80, 
011401(R) (2009a). 

[44] C.-C. Shu, K.-J. Yuan, W.-H. Hu, and S.-L. Cong, Opt. Lett. 34, 3190 
(2009b). 

[45] C. Qin, Y. Tang, Y. Wang, and B. Zhang, Phys. Rev. A 85, 053415 (2012). 

[46] S. Fleischer, Y. Zhou, R. W. Field, and K. A. Nelson, Phys. Rev. Lett. 
107, 163603 (2011). 

[47] A. G. Stepanov, L. Bonacina, S. V. ChekaUn, , and J.-P. Wolf, Opt. Lett. 
33, 2497 (2008). 



13 



